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1 Introduction 



Let X be an a-stable process X indexed by parameters (at,p), where as usual p = ¥(X\ > 0). For the 
definition and properties of stable processes and stable distributions we refer to [5], [14], [18] and [17]. 
The admissible set of parameters (a,p) is defined as 

A = {a G (0,1), p G (0,1)} U {a = 1, p = §}U {a G (1,2), p G [1 - aT 1 , a' 1 ]}. 

Note that we exclude the case when a G (0, 1) and p = 1 {p = 0}, as in this case the process X {— X} 
is a subordinator and the distribution of extrema is trivial. When a G (1,2) and p = 1 — a -1 {p = a -1 } 
the process X is spectrally positive {negative}. In this case we have complete information about the 
distribution of extrema due to the work of Bingham [6], Doney [8], Bernyk, Dalang and Peskir [4] and 
Patie [16]. 

We are interested in the distribution of the supremum of the process X, defined as 

S t = sup{X u : < u < t}. 

Note that due to the scaling property of stable processes, we have S t — t*Si, thus it is sufficient to study 
the distribution of S\. The main object of interest in this paper will be the probability density function 

p( x ) = ip(Si <x), x> 0. 
ax 

The complete description of the asymptotic behavior of p(x) as x — > + and as x — > +oo is provided 
by the following Theorem, which was proved in [13]. This Theorem generalizes results for the spectrally 
one-sided case obtained by Doney [8] and Patie [16] and the results on the two-sided case, which were 
proved by Doney and Savov [9]. 

Theorem 1. (Theorem 9, [13]) Assume that a Q. Define sequences {a m , in } m >o, n >o and {& m ,n}m>o,n>i 
as 

a = (-i) m+w f\ sin ( - {ap + 3 - 1}) n sin(7m(p + 3 - 1)} a) 

m,n r (1 - p-n- f ) T(ap + m + an) ^ sin (^) ^ sin(Traj) [) 



_ r (1 - p - n - g ) T{ap + m + an) 
r (1 + n + ^) T(-m - an) 

Then we have the following asymptotic expansions: 

P (X) ~ X ap ~ X J2Y, a ™,nX m+a \ X->0 + , (3) 

n>0 m>0 

p(x) ~ x^ a J2J2 b ^+iX~", x->+oo. (4) 



n>0 m>0 



There is a very important subclass of stable processes, for which the above result can be considerably 
strengthened. The following family of stable processes was first introduced by Doney [7]: 
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Definition 1. For k,l G Z define Ck t i as the class of stable processes with parameters (a, p) G A satisfying 

p + k = -. (5) 

a 

It turns out (see Theorem 10 in [13]), that when a process X G Ck,i for some integers k and /, then 
the coefficients a m ^ n and b m ^ n can be simplified so that they involve only finite products of lentgh not 
greater than k or /. Also, in this case the coefficients are well defined even for rational a. Moreover, 
when a G (0,1) {a G (1,2)}, the series in the right-hand side of (4) { (3) } converges to p(x) for all 
x > 0. 

It is not hard to see that parameters (a, p) G A which satisfy relation (5) form a dense set of Lebesgue 
measure zero. Therefore, Theorem 10 in [13] gives a convergent series representation for p(x) for a dense 
set of parameters. The main goal of this paper is to extend this result and to prove that in fact we have 
a convergent series representation for p(x) for almost all parameters (a,p), except when a G Q or when 
a can be approximated by rational numbers extremely well. 



2 Main Results 



The following set of real transcendental numbers, which can be approximated by rational numbers 
extremely well, was introduced in [13]. 

Definition 2. Let C be the set of all real irrational numbers x, for which there exists a constant b > 1 
such that the inequality 



p 

x 

q 



< 



l 



(6) 



is satisfied for infinitely many integers p and q. 



It is clear that £ is a proper subset of the set of Liouville numbers, which are defined by the following, 
weaker condition: for all n > 1 the inequality 



x 



< 



q r 



is satisfied for infinitely many integers p and q. A celebrated result by Liouville states that any algebraic 
number is not a Liouville number, but this is also true for many other numbers. In fact, almost every 
number is not a Liouville number, as the set of Liouville numbers, while being dense in M, has zero 
Lebesgue measure (see Theorem 32 in [12]). Therefore, the Lebesgue measure of L is also zero. 

As we will see later, the structure of the set £, can be described in terms of continued fraction 
representation of real numbers. We present here the essential results from the theory of continued 
fractions, which will be needed later. 

The continued fraction representation of a real number x is defined as 



x = [a ; at, a 2 , . . . ] = a + 



1 



a 2 
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where a® G Z and O; £ N for z > 1. For x ^ Q the continued fraction has infinitely many terms; 
truncating it after n steps gives us a rational number p n /q n = [a ; ai, a 2 , a n ], which is called the n-th 
convergent. The sequence of coprime integer numbers p n and q n can be computed recursively 



a n p n _i +p n _ 2 , p-i = 1, P-2 = 0, 
a n g n _i + g n _ 2 , <?_i = 0, g_ 2 = 1. 



(7) 



Theorem 17 in [12] tells us that the n-th convergent p n /qn provides the best rational approximation to 
x in the following sense 



Pn 




P 


x 


= min | 


x 


q n 




q 



: p G Z, 1 < q < q n \ . 



There is also a converse result (see Theorem 19 in [12]): if integers p and q satisfy 



x 



< 



2q 2 



(8) 



then p = p n and q = q n for some n. The error of the best rational approximation is bounded from above 
and below as follows 



1 



q n (q n + q n +i, 



< 



X 



Pr, 

q n 



< 



l 



qnqn+i 



(9) 



see Theorems 9 and 13 in [12]. 

The next proposition gives us an insight into the arithmetic structure of the set C We will use the 
following notation: {x} G [0, 1) denotes the fractional part of x, i.e. the distance from x to the closest 
integer not greater than x; and (x) = min({x}, {—x}) = min(|a; — n\ : n G Z), i.e. the distance from x 
to the closest integer. 

Proposition 1. 

(i) If x G C then zx G C and z + x G C for all z G Q \ {0}. 
(ii) x G C if and only if x~ l G C. 

(Hi) Let x = [do; ai, 02, ■ • •]• Then x G C if and only if there exists a constant b > 1 such that the 
inequality a n+ i > b Qn is satisfied for infinitely many n. 



(iv) x ^ £UQ if and only if 



lim 

q— >+oo 



ln(gx) 

q 



(10) 



Proof. Statement (i) follows immediately from the Definition 2. Statement (ii) also can be derived from 
the Definition 2, however it easily follows from (iii) due to the following simple property of continued 
fractions: if x > 1 and x = [ao] a\,a2, ■ ■ ■], then x~ x = [0; ao, Oi, . . .]. Thus we only need to prove (iii) 
and (iv). 

Let us prove the "if" part of (iii). Assume that there exists b > 1 such that 



a n+ i > b qn 



(11) 
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for infinitely many n. Let us consider such an index n. Using (11) and the recurrence relation (7) we 
find that q n+x > q n a n+ i > q n b qn , which together wih (9) implies that 



x 



Pn 



1 1 1 

< < -^r~ < i— 

q n q n+ i qfo q " b qn 



Since the above inequality is true for infinite many n, we conclude that x G £. 

Now, let us prove the "only if" part of (iii). Let us assume that x G £. Then there exists a constant 
b > 1 such that the inequality (6) is satisfied for infinitely many pairs (p, q). When q is large enough 
we have b~ q < l/(2g 2 ), thus p and q satisfy (8) and we conclude that the pair (p, q) must coincide with 
some continued fraction convergent (p n ,q n ). Therefore, from (9) we find 



q n (q n + q n +i) 



< 



x 



Pn 
Qn 



< 



Using the recurrence relation (7) and the above inequality we conclude that 



1 



< 



1 



1 



ql(a n+1 + 2) q n (q n + a n+1 q n + q, 



n—lj 



q n (q n + 



< 



therefore 



It is clear that if q n is large enough then 



b qn r- 
— -2>(Vb) 9 », 

q n 

Thus we have found infinitely many indices n such that that a n+ i > (Vb) qn . This ends the proof of (iii). 

Let us prove the "if" part of (iv). Since i^£UQwe know that for every b > 1, all q G N sufficiently 
large and all p G Z we have 



\qx 



i 9 1 

Fl ¥ b q 



therefore ln(gx) > —qhx(b) for all q sufficiently large, which shows that 

liminf^M>0. 



q— >+oo 



q 



Since we also have \n(qx)/q < for all q, it implies condition (10). The "only if" part of (iv) can be 
verified in exactly the same way and we leave the details to the reader. □ 



Proposition 1 shows that the set C has quite an interesting structure. First of all, property (iii) gives 
us a simple method to construct a number x G £ (just define recursively a n+ \ = 2 Qn ), therefore this set 
is not empty. The set £ is closed under addition and multiplication by rational numbers, therefore it is 
dense in R. It is a subset of the set of Liouville numbers, therefore it has Lebesgue measure zero. In 
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fact, the Hausdorff dimension of the set of Liouville numbers is also zero (see Theorem 2.4 in [15]), and 
of course the same is true for L. 

We would like to stress that the elements of the set L are rather unusual transcendental numbers with 
quite extreme arithmetic properties. The reason for this lies in propety (iii) of the above proposition. 
Assume that x is an irrational number. Theorem 12 in [12] tells us that for every irrational number x 

<?n>2^, n>2, 

therefore using Proposition 1, (iii) we find that for every number x G C there exists a constant b > 1 
such that the inequality 



a n > 



(12) 



is satisfied for infinitely many n. Yet one can prove that for all x except a set of Lebesgue measure zero 
and all b > 1 



a n = 0(b n ), n -> +oo. 



(13) 



This result follows from the famous Levy-Khinchin theorem, which states that for almost all x it is true 
that ^fq^ — > 7 as n — > +oo, where 7 = exp(-7r 2 /(12 ln(2))) (see Theorem 31 and the footnote on page 66 
in [12]). We see now that there is a very large gap between the growth rate of coefficients a n for elements 
of C given by (12) and almost all other real numbers given by (13). This shows that the elements of the 
set L should be rather exceptional numbers. 

Our main result in this paper is the following Theorem, which gives an absolutely convergent series 
representation for p(x) for every irrational a which is not in the set C 



Theorem 2. Assume that a ^ C U Q. Then for all x > 



p(x) 



x 



1 a $Z5Z &m - nH 



n>0 m>0 



X 



OLp- 



rn+an 



if a e (0,1), 
if a G (1,2), 



(14) 



n>0 m>0 

where a m ^ n and b m ^ n are defined by (1) and (2). 

Before we prove Theorem 2, let us establish the following result. 
Lemma 1. Assume that x ^ C U Q. Then 

k 



J|sec(vr/x)| = 2 fc+o(fc) , k ->■ +00, 
1=1 

k 

J|csc(7r/x)| = 2 fc+o(fc) , k -> +00. 



(15) 
(16) 



1=1 



Proof. Let us prove (15). We use the following easily verified fact 



{y}-- 2 



1 \ <2y> 
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and Proposition 1 (iv) to check that for all x ^ C U Q 



Urn H{**}-a =a 

Since for all x ^ Z + | it is true that 

, M 2 



Sec 7TX < 



we conclude that 



In I sec(7rA^x)| 
lim — ! £ ^ = 0. 

Applying Theorem 1 from [3] we see that the above statement implies 

k 



lim — y^ln I sec({7r/x})| = ln(2), 



fc->+oo k 

1=1 



which is equivalent to (15). Formula (16) can be verified in the same way; it also follows from Theorem 
2 in [11]. □ 



Using the results of Lemma 1 it is not hard to prove that both series in the right-hand side of (14) 
converge for all x > 0. However this does not guarantee that the series converges to p(x). As an example, 
consider the asymptotic expansion 

f(x) = e~ x -\ — x -> +oo. 

v ' x -1 ^ 

n>l 

The asymptotic series in the right-hand side of the above equation converges for all x > 1, but the limit 
is not equal to f(x). Therefore, in order to prove Theorem 2 we would need to do some more work. 
Our main tool will be the Mellin transform of Si, which is defined as 

M(w) = E[S7 _1 ] = J p(x)x w ~ 1 dx, Re{w) = 1. 

This function was studied in [13], where it was proved that it can be analytically continued to a mero- 
morphic function. There exists an explicit expression (Theorem 8 in [13]) for A4(s) in terms of the 
double gamma function, see [1] and [2] for the definition and properties of the double gamma function. 
It is also known that A4(s) satisfies several functional equations (see Theorem 7 in [13]) and if a £ Q 
and X £ Ck,i it has simple poles at the points 

{ S m,n\m>l,n>l = {m + an} m >l >n >l, {s~ ;n } m > ,n>0 = {1 ~ Oip - TTl - Cm} m > , n > , 

with residues given by 

Res (M(s) : s+ „) = -b m -i, n , Res (M(s) : s~ n ) = a m , n , 
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see Lemma 2 in [13]. 



Proof of Theorem 2: In what follows we will always use the principal branch of the logarithm, which 
is defined in the domain | arg(z)| < ir by requiring that ln(l) = 0. Similarly, the power function will be 
defined as z a = exp(aln(z)) in the domain | arg(z)| < n. 

Assume that a ^ £ U Q, a £ (1, 2) and define = 1 — ap + f — k. Note that c\. ^ s~ n for all 
k,m,n > 0; if this was not true then a would be a rational number. Therefore Ai(s) does not have 
singularities on the lines + iR for all k > 0. 

Following the proof of Theorem 9 in [13], we start with the inverse Mellin transform representation 



1 
2vri 



M(s)x s ds, x > 0. 



The integral in the right-hand side of the above formula converges absolutely, since Mis) — > exponen- 
tially fast as Im(s) — > oo, see Lemma 3 in [13]. Shifting the contour of integration 1 + iR i— >■ + iR and 
taking into account the residues at the poles s = s~ , we find that 



p(x) = ^Res(A^(s) : s~ 



x x 



.,n _|_ 



2m 



M(s)x Ms. 



(17) 



where the summation is over all m > 0,n > 0, such that s~ n > Now our goal is to prove that as 
k — > +oo, the integral in the right-hand side of (17) - which we will denote by Ik{x) - converges to zero 
for every x > 0. 

First, we perform a change of variables 

I k (x) = \ J M(s)x~ s ds = x ~ 1+ap ~^ +k J M (l - ap + | - k + iu) x~ iu du. (18) 

c fe +iR R 

Using equation (6.2) in [13] and the reflection formula for the Gamma function (see formula (8.334.3) in 
[10]) we find that 



Mis) 



Mis 



r( s + i)r 



!-('+!) 



sing(l- S )) 



sin(^(ap 



0)' 



seC. 



Iterating the above identity k times gives us 



M(s] 



M(s + k) 



r( 8 + k)r (i=*=2) 



n 



sin (^(2- j -s)) 
Lsin(^(ap-2+j + S ))' 



sGC. 



(19) 



Next, we use the fact that ap < 1 for («,p) £ ^4, therefore —1 + ap — ^ < while for k > 2 we have 
— 1 + ap — 7j + k > 0, which in turn implies the estimate 

x -l+ap-f +k < ^ + ^ fc > 2j ^ > q ^ 

We express M(s) in terms of .M(s + k) with the help of (19), substitute the resulting expression into 
(18) and use (20) to obtain the following estimate 



\h(x)\ < 



[l + xf J \M (l-ap+ — + iiij | x |Fi(tt;fc)| x \F 2 (u; k)\du, k > 2, x > 0, 



(21) 



where we have denoted 



Fx («;*)=] 



sin(^ (ap-2-iu + k+l-j)) 



(22) 



and 



F 2 (u;fc) 



r(l-ap+f + m - fc)r (p - 1 - ig + g) 

r(i- ap + f + iu )r(p-i-i2) 



(23) 



First, let us estimate Fi(u;k) as fc — > +00. Our goal is to prove that as — >• +00 the function 
Fi(u; k) does not grow faster than some exponential function of k, and that this bound is uniform in 
!t£l. From the trigonometric identities 

sm(x + iy) = sin(x) cosh(y) + i cos(x) sinh(y), 
I sin(x + iy)\ 2 = cosh(y) 2 — cos(x) 2 

we see that | sin(x)| cosh(y) < | sin(x + iy)\ < cosh(y), therefore 



sin(a + iy) 



< 



sin(6 + iy) 

Combining the above estimate with (22) we see that 

k 



sin(6)| 



\Fi(u]k)\ < Y[ 



1 71 f a 



n 

1=1 



sec 



(24) 



where in the last step we have changed the index of summation I — k + 1 — j. As we have established 
in the Proposition 1, a ^ C U Q implies o _1 ^ £UQ. Therefore the inequality (24) and Lemma 1 tell 
us that there exists a constant C\ = C\{a) > such that 



< Ci3 fc , k>2, MGR. 



(25) 



Now we will deal with F 2 (u; k) defined by (23). Our goal is to prove that as k — > +00 this function 
converges to zero faster than any exponential function of k, and that this happens uniformly in u e M. 
More precisely, we will prove that there exists a constant C 2 = 6*2(0;) > such that 



\F 2 (u; k)\ < C 2 (l + M)exp 



a 



a 



■(k-l)ln(k-l)+2k\A\ , k > 2, ue 



(26) 



where A = (1 + ln(a) — a) /a. 

Using the reflection formula for the gamma function we rewrite F 2 (u; k) as 



F 2 ( U] k) = (-l) k g(u)x 



r (p-i-^+a 

r [ap — f — iu + k) 



where 



9W 



71 



a 



sin ( 7r ( ap — — — iu ) ) T ( 1 



a 



ap + — + iu ) r I p — - — i— 



u 



o 
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Next we use the following asymptotic expression (see formula (8.328.1) in [10]) 

\T(x + iy)\ ~ \Fln\y\ x ~^e~ 1 ^ y \ y ->• oo, 
and the fact that | sin(x + iy) \ ~ exp(|y|)/2 as y — > oo to obtain 



\9W 



7T 



2 v i i 



p-i 



n -1 



_ ZL M 

e 2 a 



(27) 



P-1L.|(«-1)M) P -( Q -D^ 



ar \u 



as m — 7- oo. We use the inequality (a — 1) (p — |) < a — 1 < 1 and the fact that g(w) is continuous in w 
to conclude that there exists a constant C3 = C^ia) > such that 



<?(w) < C 3 (l + |M|)e- (a_1) ^ 



u e 



(21 



Next, from the Stirling's asymptotic formula for the Gamma function (see formula (8.237) in [10]) we 
find that 



rls) 



a cxp 



s — y 00, Re(s) > 0, 



where A = (1 + ln(ct) — a) /a (note that A < 0). Again, the function in the left-hand side of the above 
equation is continuous in the half-plane Re(s) > 0; this fact and the above asymptotic formula imply 
that there exists a constant C4 = CAa) > such that 



r(S) 



r( a ) 



< c 4 



exp 



I - — -ln(s) + A 



Re(s) > 0. 



We set s = ap — ^ — iu + k (note that for k > 2 we have Re(s) > 0) and use (28) and the above estimate 
to conclude that for k > 2 



\F 2 {u-k)\ = \g{u)\x 



r( s ) 

< C 5 (l + |n|)exp 
= C 5 (l + |u|) exp 

< C 5 (l + |u|) exp 



(29) 



1)5^ + Re <! -s ( ln(s) + A 



fa-T 



2a 

7l\u\ 

~2a 



a 



- Re(s) ( - — - h 



a 



a-1 

u arg(s) 

a 



-Re(s) ( - — - h 



a 



where C5 = C3 x C4 and in the last step we have used the fact that | arg(s)| < tv/2. Using the following 
facts: (i) \s\ > Re(s) > k — 1 and (ii) Re(s) < 2k we check that 



cxp 



, a — 1 
-Re(s) ( In \s\ + A 



a 



< exp 



a — 1 



a 



{k - 1) In(fc - 1) + 2k\A\ 



(30) 



and combining (29) and (30) we finally conclude that (26) is true with C2 = C$ 
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Combining (21), (25) and (26) we see that 

\h(x)\ < C 6 (3(l + x)e 2 ^) fc e-^ (fc - 1)ln(fe - 1) , k > 2, 

where 

C 6 = C 5 x Ci x J \M (l - ap + ^ + iuj | x (1 + |u|)du. 

R 

Therefore Ik{x) — > as k — > +oo and we have proved the second series representation in (14). The 
proof of the first series representation in (14) is identical, except that we have to shift the contour of 
integration in the opposite direction. 

□ 

Theorem 2 gives us a convergent series representation for p(x) provided that a £ £UQ. Let us 
consider what happens in the remaining cases. When a G Q and X ^ Ck,i for all then Theorems 

1 and 2 can not be valid for the following reason. Using Theorem 8 in [13] and properties of the double 
gamma function (see [1] and [2]) we conclude that in this case the function Ai(s) will have multiple poles 
on the real line. This means that the asymptotic expansions (3), (4) would include terms of the form 

Res{M(s) : s = w) x ln(x) i x""', < j < k - 1, 

where w is the pole of A4(s) of multiplicity k. This also means that the coefficients of such an expansion, 
given by the residues of Ai(s), will depend on the derivatives of the double gamma function, and in this 
case it will be much harder (maybe even impossible) to evaluate them analytically. 

In the remaining case, when a 6 C, it is not clear what happens with the series representation for 
p(x) given in (14). By inspecting the form of coefficients a mj „ and h m , n defined in (1) and (2) we see that 
the series in (14) might still converge, as the growth of the product of sin(-) functions in the numerator 
in (1) might "cancel" the growth of the denominator. To prove the convergence one would need some 
results on the asymptotic behavior as k — > +oo of the following product 

k 

Y[ I csc(tt(/x + y))\. 
i=i 

Results in [3] give us information about these products in the case when y G Q and x ^ C, but we were 
not able to find any general results in the existing literature. In any case, even if one were able to prove 
convergence of the series (14) for a more general class of parameter a, it would not imply that the series 
converges to p(x), see the discussion on page 7. To summarize, when a G C the question of validity of 
Theorem 2 is open, and there are three possible scenarios: (i) the series in (14) do not converge, (ii) they 
converge, but not to p(x), (iii) they converge to p(x). 
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